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1 Introduction 



In minimal supersymmetric extensions of the Standard Model (SM), soft Supersymmetry (SUSY) 
breaking terms [1] are introduced to break SUSY without spoiling the cancellation of quadratic 
divergences in the process of renormalization. These terms must have dimensionful couplings, 
whose values determine the scale M SUSY , lower than a few TeV, above which SUSY is restored; 
they are also responsible for the mass splittings inside the supersymmetric multiplets. Little is 
known for sure about the origin of these SUSY-breaking terms. Upcoming accelerators will test 
the energy range where we hope that the first supersymmetric particles will be found. From their 
masses and couplings we could learn about the pattern of SUSY-breaking at low energies, which 
translates, through the renormalization group equations, into the pattern of breaking at the scale 
at which SUSY-breaking is transmitted to the observable sector. Meanwhile, one can obtain some 
information on the soft terms by looking at any low-energy observables sensitive to their values, 
and in particular to the Yukawa sector of the theory. 

In this work we consider the simplest supersymmetric version of the SM, the Minimal Super- 
symmetric Standard Model (MSSM) [2]. We analyse the limit of a large ratio V2/V1 = tan/3 of the 
vacuum expectation values v\, V2 of the Higgs doublets. We show that in this limit a large class 
of physical observables involving the Yukawa coupling of the physical charged Higgs boson can be 
described in terms of a two-Higgs-doublets model (2HDM) [3] effective lagrangian, with specific 
constraints from the underlying MSSM dynamics. 

The finding of a charged Higgs boson would be instant evidence for physics beyond the SM. 
It would also be consistent with low-energy SUSY, as all supersymmetric extensions of the SM 
contain at least a charged Higgs boson, H^. Current experiments, looking at the kinematical 
region M H+ < m t — m b , have been able to place an absolute bound of M H+ > 71.0 GeV at the 
95% confidence level [4] and/or to exclude regions of the M H+ -BlZ{t — > bH + ) plane [5, 6]. 1 If 
the charged Higgs mass happens to be greater than the top mass, future e + e~, pp and even e~p 
accelerators will have a chance to find it [8-10]. 

Present bounds from LEP on a SM light Higgs boson, M h sm > 105.6 GeV [11], are beginning 
to put strong constraints on values of tan/3 lower than a few, a region that can only be consistent 
with low-energy SUSY if the third-generation squark masses are large, of the order of a TeV and, 
in addition, if the mixing parameters in the stop sector are of the order of, or larger than, the 
stop masses. Therefore, the LEP limits give a strong motivation for the study of the large tan/3 
region. The region of large values of tan/3 is also theoretically appealing, since it is consistent 
with the approximate unification of the top and bottom Yukawa couplings at high energies, as 
happens in minimal SO(10) models [12, 13]. The aim of this work is to compare, for large values 
of the tan/3 parameter, the effective potential results truncated at one loop with the diagrammatic 
one-loop computation for the supersymmetric QCD (SUSY-QCD) and electroweak (SUSY-EW) 
corrections in the coupling of ibH + [14-17]. We then use the effective potential approach to include 
a resummation of the SUSY-QCD and SUSY-EW effects and we show how relevant these higher- 
order effects are to the final evaluation of the H + — ► tb and t — > b H + partial decay rates. 

Although diagrammatic computations of the 0(ga s ) quantum corrections to these observables 
have existed in the literature for several years, either in the context of a generic two-Higgs-doublets 
model [18-21], 2 or in supersymmetric extensions of the Standard Model [14-17], our analysis goes 
beyond these studies in the following: 

• It resums leading and next-to-leading logarithms of the type a s log(m b /m t ) or 

1 See also the study in ref. [7], where it is shown how these bounds are affected by some usually overlooked decay 
modes in the intermediate tan/3 > 1 region. 

2 For the QCD corrections to the neutral Higgs decay rate the reader is referred to [22, 23] and references therein. 
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Figure 1: One- loop SUSY-QCD diagram contributing to the effective coupling Ah b . The solid lines 
inside the loop denote the gluino propagator, the dashed lines correspond to sbottom propagators. 
The cross represents the Mg insertion coming from the gluino propagator. 

a s log(m b /M H+ ), because these terms are of the same size as the tree-level result. 

• It includes the potentially large supersymmetric corrections responsible for the leading be- 
haviour at large tan/3 > 10 values, with an improved treatment of the higher-order contribu- 
tions incorporated into the effective lagrangian: the corrections of order (a s fi tan (3 /M SUSY ) n 
are included to all orders n = 1, 2, . . . These corrections do not vanish if the parameter fi and 
the soft SUSY-breaking masses are pushed to large values, which is a reflection of the lack of 
supersymmetry in the low-energy theory. 

• It is well suited for numerical evaluation, because it includes all the relevant terms, by means 
of very simple formulae. Therefore, the bulk of the quantum corrections can be implemented 
in a fast Monte Carlo generator. 

We would like to stress the second point: even for a heavy supersymmetric spectrum, depending 
on the ratios and relative signs of the Higgs mass parameter, fj,, and of the soft SUSY-breaking 
parameters involved, the supersymmetric QCD and EW corrections can be very large, a situation 
in which the higher-order effects are sizeable. 

The text is organized as follows. In section 2 we derive the coefficients of the 2HDM effective 
lagrangian which are affected by large SUSY threshold effects. Section 3 provides simple analytical 
expressions for the QCD and electroweak quantum corrections to the t b H + and H + — ► 1 6 partial 
decay rates, including the resummation of the large leading and next-to-leading QCD logarithms 
and of the potentially large tan/3-enhanced SUSY corrections. Section 4 is devoted to the numerical 
analysis of the partial widths, comparing them to the previously existing one-loop results [15, 17]. 
To exemplify the importance of these novel computations, we show in section 5 their effects on the 
BWs of t — ► b H + and H + — ► tb. As an example we study the effects of these results on the limits 
on the M H+ mass derived by the DO collaboration (similar limits have been obtained by the CDF 
collaboration) via the indirect search of the charged Higgs in tt decays. We reserve section 6 for 
our summary and conclusions. 

2 Effective lagrangian 

2.1 Supersymmetric corrections 

The effective 2HDM lagrangian contains the following couplings of the bottom quark to the CP-even 
neutral Higgs bosons [24]: 

h b H°bb + Ah b H$bb . (1) 
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The H^bb tree-level coupling is forbidden in the MSSM. Yet a non-vanishing Ah b is dynamically 
generated at the one-loop level by the diagram of fig. 1. 

Although Ah b is loop-suppressed, once the Higgs fields Hf 2 acquire their vacuum expectation 
values vip, the small Ahb shift induces a potentially large modification of the tree-level relation 
between the bottom mass and its Yukawa coupling, because it is enhanced by tan/3 = V2/V1: 

m b = hbV\ — ► m b = v\ (h b + Ahb tan/3) = h^vi (1 + Am b ) . (2) 

Since the numerical value of mb is fixed from experiment, equation (2) induces a change in the 
effective Yukawa coupling. This affects not only the CP-even neutral Higgs field, but the whole 
Higgs multiplet, with phenomenological consequences for the charged Higgs particle. In particular, 
eq. (2) modifies the Yukawa coupling of the charged Higgs to top and bottom quarks as follows: 

h b sm(3 = ^ tan/3 — ► h b = — tan/3, (3) 

v v 1 + Am b 



where v = \Jv\ + v\ ~ 174 GeV. In the last equation we have assumed a large tan/3 regime. 

It turns out that, in the MSSM with large tan/3, the dominant supersymmetric radiative cor- 
rections to the Yukawa interactions of the Higgs doublet H\ = (H± \H®) stem from the relation 
(3). Explicit loop corrections to the H\ff Yukawa coupling are suppressed by at least one power 
of tan/3. This remarkable feature has far-reaching consequences: first in observables involving the 
coupling hb of H\ to bottom quarks the MSSM behaves like a two-Higgs-doublets model. The main 
effect of a heavy SUSY spectrum is to modify the coupling strength via Am b in eq. (2), which 
depends on the masses of the supersymmetric particles. In certain regions of the parameter space 
a sizeable enhancement of h b occurs. Secondly these dominant corrections encoded in Am b are 
universal. They are not only equal for the neutral and the charged Higgs bosons, on which we will 
focus in the following, but they are also independent of the kinematical configuration. This means 
that they affect the decay rate of a charged Higgs into a top and bottom (anti-) quark in the same 
way as the tbH + vertex in a rare 6-decay amplitude or, after replacing the top by a charm quark, as 
Higgs-mediated b — > c decays. Further the universality property of these tan/3-enhanced radiative 
corrections allows for a simple inclusion into the Higgs search analysis. 

The proper tool to describe such universal effects is an effective lagrangian. Expanding (1) to 
include the charged Higgs sector one finds that the relevant terms in the large tan/3 limit are: 

C = -h b b L b R H^ + h b V tb sm(3t L b R H + - Ah b b L b R H% + h.c. (4) 

Ahb is the loop-induced Yukawa coupling associated with the supersymmetric QCD corrections in 
fig. 1 and similar electroweak contributions. H + is the physical charged Higgs boson. The Higgs 
mechanism defines the relation between the bottom mass m b and the couplings h b and Ahb m 
C: calculating the tree-level bbHf and one-loop bbH® vertices with zero Higgs momentum, and 
replacing the Higgs fields by their vacuum expectation values v± t 2, yields the desired relation in 
eqs. (2) and (3): 



Amb = ^ tan/3 = Am S b QCD + Am s b EW , (5) 
ti b 



which contains the tan/3-enhanced radiative corrections. The supersymmetric QCD corrections of 
fig. 1 read [13] 

Am SQCD = ^± M ^ tan p I{m ^ m ^ Mj .^ (6) 



3 There are similar diagrams involving supersymmetric electroweak quantum corrections, see section 3.2. 
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Here a s is the strong coupling constant and fi is the mass parameter coefficient of the eijH\H 2 
term in the superpotential. The vertex function /, which depends on the masses m~ of the two 

01,2 

bottom squark mass eigenstates and the gluino mass Mg, reads [13] 

/( °' 6 ' C) = (a 2 -b 2 )(b 2 -c 2 )(a 2 -c 2 ) (f# log + #<? h ~ 2 + c 2 a 2 log ^ . ( 7) 

An interesting limit of eq. (6) applies when all mass parameters are of equal size. One has, depending 
on the sign of fi 

Am S b QCD = ± a s{Q = M SUSY ) tan/3 ^ 

clearly showing that the effect does not vanish for a heavy SUSY spectrum and can be of 0(1) for 
large tan/3 values. 

For sizeable values of the trilinear soft SUSY-breaking parameter A t , the supersymmetric elec- 
troweak corrections are dominated by the charged higgsino-stop contribution, which is proportional 
to the square of the top Yukawa coupling, ht = m t /v 2 . Wino-sbottom contributions are generally 
smaller, being proportional to the square of the SU(2)l gauge coupling, g, and to the soft SUSY 
breaking mass parameter M 2 . Neglecting the bino effects, which we found to be numerically irrel- 
evant, these corrections read [25] 

h 2 

Amf EW = j^nA t tanpI( mii ,m h ,n) 



Q 2 

/ uM 2 tan/3 



16vr 2 



cos 2 9 i i"(m f ,M 2 ,fJb) + sin 2 % i"(m ? ,M 2 ,fi) 



+ ^co S 2 9 i I(m ii ,M 2 ,^ + ^ S m 2 9~ b I(m- b2 ,M 2 ,fi)] . (9) 

When including radiative corrections, one has to specify the definition of the quark mass m b 
appearing in the leading order: m b denotes the pole mass corresponding to the on-shell renormal- 
ization scheme, in which the on-shell self-energy is exactly cancelled by the mass counterterm. 

Note that the supersymmetric corrections contained in Am b enter hb in eq. (3) as a factor 1/(1 + 
Am 6 ). To order a s one is entitled to expand this factor as (1 — AmJ. In the phenomenologically 
most interesting case of a large |AmJ of 0(1), this leads to disturbingly large numerical ambiguities. 
Their resolution seems to require painful higher-order loop calculations, and a large |Am 6 | may 
even put perturbation theory into doubt. Yet these tan/3-enhanced contributions have the surprising 
feature that they are absent in higher orders: 

There are no contributions to Am b of order 

{<XsJT—^P\ (10) 

\ 1V1 SUSY J 

for n > 2. 

Here M SUSY represents a generic mass of the supersymmetric particles. An analogous result 
applies to the electroweak corrections. In other words, to the considered order, Am b is a one-loop 
exact quantity, and the factor 1/(1 + Am b ) contains the corrections to hb of the form in (10) to all 
orders in a s . 

To prove our theorem, consider possible n-loop SUSY-QCD contributions to Am b proportional 
to tan n (3: the only possible source of additional factors of tan/3 is the off-diagonal element of 
the bottom squark mass matrix, —/j,m b tanf3, which can enter the result via the squark masses 
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y Figure 3: One-loop diagram derived from the 

effective lagrangian in (4) corresponding to 
the diagram in fig. 2. It contains the large 
Figure 2: Two-loop SUSY-QCD di- logarithm of fig. 2 as \ogm b /Q. This loga- 
agram containing a large logarithm rithm is summed to all orders by solving the 

log (M SUSY /m b ). renormalization group for hi in eq. (4). 

as m- 2 — m- 2 ~ 2m, a tan/3 or through counterterms to the squark masses. It is easier to track 

02 01 u 

the factors of fim b tan/3 by working with "chiral" squark eigenstates and assigning these factors to 
"chirality flipping" two-squark vertices. Thus any extra factor of tan/3 is necessarily accompanied 
by a factor of m b fi. This dimensionful factor is multiplied with some power of inverse masses 
stemming from the loop integrals. The next step in our reasoning is to show that the loop integrals 
always give powers of 1/M SUSY and can never produce a factor of l/m b . The appearance of any 
inverse power of m b in a loop integral would imply a power-like infrared singularity in the limit 
m — > with gluino and squark masses held fixed. But the KLN theorem [26] guarantees the 
absence of any infrared divergence in all bare diagrams except for those in which gluons couple to 
the 6-quark lines. A two-loop example of the latter set is shown in fig. 2. The infrared behaviour 
of these diagrams can be studied with the help of the operator product expansion (OPE). The 
result of the OPE is nothing but the effective lagrangian in (4). To apply the OPE to our problem 
we first have to contract the lines with heavy super symmetric particles to a point, i.e. we replace 
the MSSM by an effective theory in which the heavy SUSY particles are integrated out. For the 
case of the diagram in fig. 2 this yields the diagram in fig. 3, in which the loop- induced interaction 
is represented by the dimension-4 operator bbH®. The information on the heavy SUSY masses is 
contained in the Wilson coefficient A/i in eq. (4). The key feature of the OPE exploited in our 
proof is the fact that the effective diagram in fig. 3 and the original diagram in fig. 2 have the 
same infrared behaviour. Power counting shows that the diagram of fig. 3 has dimension zero. It 
depends only on m b and the renormalization scale Q. Since Q enters the result logarithmically, 
the diagram of fig. 3 depends on m b as logm b /Q, no power- like dependence on m b is possible. 
This argument — essentially power counting — immediately extends to higher orders. Terms from 
diagrams in which gluons are connected with the 6-quark line and one of the SUSY-particle lines 
in the heavy loop, are either infrared-finite or suppressed by even one more power of m b /M SUSY , 
because they are represented in the OPE by operators with dimension higher than 4. Finally there 
are diagrams with counterterms. In mass-independent renormalization schemes the counterterms 
are polynomial in m b . In the on-shell scheme the diagrams with counterterms can be infrared- 
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Figure 4: One-loop SUSY-QCD diagram contributing to Am b . 

divergent for m b — > 0, but only logarithmically. In conclusion the loop integrals cannot give factors 
of l/m b . Therefore any correction to Amj of order a™tan n /3 comes with a suppression factor of 
m b /Mg USY - Higher-order loop corrections to Am b are therefore either suppressed by m b /M SUSY 
or lack the enhancement factor of tan/3, which proves our theorem. 

So far we have discussed Am b from the one-loop vertex function of fig. 1 as in [13]. A different 
viewpoint has been taken e.g. in [14]: the renormalization of the Yukawa coupling (m b /u) tan/3 is 
performed by adding the mass counterterm to m b . In the large tan/3 limit and to one-loop order, 
this amounts to the replacement 

^ tan/3 — ► h b = ^ (1 - Am h ) tan/3 (11) 
v v 

instead of (3). This procedure gives the correct renormalization of the Yukawa coupling in regular- 
ization schemes respecting gauge symmetry [27], such as dimensional regularization. The relation 
to the Yukawa renormalization using the vertex function in fig. 1 leading to (2) is provided by a 
Slavov- Taylor identity [27]. In general a correction factor related to the anomalous dimension of 
the quark mass occurs in (2), but the large tan/3-enhanced contributions considered by us are finite 
and do not contribute to the anomalous mass dimension. To one-loop order, eqs. (11) and (3) 
are equivalent. Yet the crucial difference here is the point that — Am b in eq. (11) stems from the 
supersymmetric contribution to the quark self-energy diagram in fig. 4. While the vertex diagram 
has dimension zero, the self-energy diagram has dimension one and the above proof does not apply. 
Indeed, higher-order corrections to fig. 4 do contain corrections of the type in (10). In Appendix 
A these corrections are identified and it is shown that they sum to 

(12) 



1 + Am b ' 

so that both approaches lead to the same result (3) to all orders in (n/M SUSY ) a s tan/3. 
2.2 Renormalization group improvement 

The tan/3-enhanced supersymmetric corrections discussed so far are not the only universal correc- 
tions. It is well known that standard QCD corrections to transitions involving Yukawa couplings 
contain logarithms log(Q/m b ), where Q is the characteristic energy scale of the process. For the 
decays discussed in sects. 3-5 one has Q = m H + or Q = m t and a s log (Q/m b ) is of 0(1) thereby 
spoiling ordinary perturbation theory. The summation of the leading logarithms 

<log n -^-, n = 0,l,2... (13) 

m b 

to all orders in perturbation theory has been performed in [22] for the standard QCD corrections 
to the tibRH + Yukawa interaction. This summation is effectively performed by evaluating the 
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running Yukawa coupling hb at the renormalization scale Q. This amounts to the use of the 
running mass at the scale Q, m b (Q), after expressing h b sin f3 in terms of (m b /v) tan/3. Hence these 
large logarithms are likewise universal, depending only on the energy scale Q at which the Yukawa 
coupling is probed, and can also be absorbed into the effective lagrangian. 

The full one-loop QCD corrections to neutral [22] and charged [18] Higgs decay and top decay 
[21] also contain non- logarithmic terms of the order a s . A consistent use of these one- loop corrected 
expressions therefore requires the summation of the next-todeading logarithms 

a n+i log n_Q_ n = 0,1,2... (14) 

to all orders, because all these terms have the same size as the onedoop finite terms. Since squarks 
and gluinos are heavy, leading logarithms of the type in (13) are absent in the supersymmetric 
corrections shown in fig. 1. It is important to note, however, that this is no longer true for the 
next-todeading logarithms: dressing fig. 1 with n gluons leads to diagrams involving the logarithmic 
terms of (14). A two-loop diagram yielding a term of order a 2 s \og(Q/m b ) is shown in fig. 2. These 
subleading logarithms have escaped attention so far. In the remainder of this section we will address 
their proper summation. 

In [22] it has been proved that all leading logarithms occurring in neutral Higgs decays can be 
absorbed into the running mass m b (Q). This proof is based on the KLN theorem [26] and exploits 
the fact that there are only two mass scales, m b and Q, in the loop corrections to neutral Higgs 
decay. This reasoning cannot be extended to the next-to-leading logarithms accompanying the 
supersymmetric corrections of fig. 2, where both heavy and light masses occur in the loops. Here 
we will use the OPE instead and apply standard renormalization group methods to the effective 
coupling in eq. (4). This is not only much more elegant than the method used in [22], it will also 
show us how to consistently combine the summation of large logarithms with the all-order result 
of the tan/3-enhanced terms derived in section 2.1. 

To apply the OPE and the renormalization group one must first employ a mass-independent 
renormalization scheme, such as the MS scheme [28]. At the scale Q = M SUSY the heavy particles, 
squarks and gluinos, are integrated out. The interaction mediated by the loop diagram in fig. 1 is 
now represented by the effective operator bbH®. Its Wilson coefficient equals 

-Ah b (Q = M SUSY ). (15) 

Here and in the following, MS quantities are overlined. The renormalization scale Q is explicitly 
displayed in (15). Note that Ahb depends on Q through a s , M~ g and the squark masses. The relation 
(3) between hj, and m b is defined at the low scale Q = m b . Hence we must evolve (15) down to 
Q = m b . Since we encounter the same operator bb as in the leading order, the renormalization 
group evolution down to Q = m b is also identical to the leading-order evolution and just amounts 
to the use of the running Yukawa coupling h b (Q = m b ) in the desired relation: 

h h (Q=m b ) = m ^Q = m b) tan/? (16) 

v 1 + Am b (Q = M SUSY ) 

Notice that Am b is evaluated at the high scale Q = M SUSY : the heavy particles 'freeze out' at 
the heavy scale Q = M SUSY and the strong coupling a s in Am b likewise enters the result at this 
scale. This can be intuitively understood, as the loop momenta in fig. 1 probe the strong coupling 
at typical scales of order M SUSY . Further any renormalization group running below Q = M SUSY 
is done with the standard model result for /3-functions and anomalous dimensions. Since the QCD 
contributions to the anomalous dimensions of hf, and are the same, h b at an arbitrary scale Q 
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is given by 

h(Q) 



m b (Q) 



v 1 + Am b (M SUSY ) 



tan/3. 



(17) 



If one expands h b (M SUSY ) around h b (m b ) to order a 2 s , one reproduces the large logarithm of the 
form log(M SUSY /m b ) contained in the diagram of fig. 2. The running mass must be evaluated with 



the next-to-leading order formula: 

rn b (Q) = U 6 (Q,m t ) ■ U 5 (m t ,m b ) ■m b (m b ) , 



(18) 



where we have assumed that there are no other coloured particles with masses between Q and m t . 
The evolution factor Uj reads 



U f (Q 2 ,Qi) = 
d (f) = 

J(» = 



1 + 



a s (Q 1 ) 
12 

33-2/' 

8982 - 504/ + 40/ 2 
3(33-2/) 2 



a s (Qi) - a s (Q 2 ) 
4tt 



jU) 



(19) 



Here / is the number of active quark flavours. For Q < mt one must replace Uq (Q,mt) - Us (mt,m b ) 
by U^{Q,m b ) in eq. (18). depends on the renormalization scheme, the result in eq. (19) is 

specific to the MS scheme. The 6-quark mass in this scheme is accurately known from T(IS') 
spectroscopy and momenta of the bb production cross section [29]: 



m b (m b ) 



(4.25 ± 0.08) GeV. 



(20) 



Physical observables such as the H + and top decay rates discussed in sections 3-5 are scheme 
independent to the calculated order. Passing to a different renormalization scheme would change 
JW, but in eq. (18) the change in a s (771,5) is compensated by a corresponding change in the 
numerical value of m b (m b ). Likewise the scheme dependence in a s (Q) is compensated by the 
one-loop standard QCD corrections [18, 21] to the decay rates. This concludes the discussion of 
the universal renormalization group effects. A discussion of additional aspects specific to the decay 
rates F(t — ► bH + ) and T(H + — > tb) can be found in Appendix B. 
Finally we arrive at the desired effective lagrangian for large tan/3: 



9 rn b {Q) 
2M W 1 + Am 6 



tanP A ibj 5 b(Q) 
+ V2V tb tan/3 H + t L b R {Q) + V2V t * b tan/3 H-~b R t L (Q) 



+ 



sin a 

cos/3 
cos a 
cos /3 



cos a , r 
Am t — ) hbb(Q) 



sin/3 



+ Am b ^)Hbb(Q) 



sin/3 



(21) 



where the renormalization scale Q entering m b and the renormalization constants of the quark 
bilinears are explicitly shown. In equation (21) we have expressed C in terms of the physical Higgs 
fields H, h, A and H + and traded v for the W mass and the SU(2) gauge coupling g. We have used 
the standard convention [3, 24] for these fields and the h-H mixing angle a. For completeness also 
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the coupling of the CP-odd Higgs boson A has been included. The phenomenology of the neutral 
Higgs bosons in the large tan/3 regime has been studied in detail in [24]. 

The effective lagrangian in eq. (21) describes the Abj^b and H + t]J)ji interactions correctly for 
large tan/3, irrespective of the mass hierarchy between M SUSY and M H+ . Even if M SUSY ~ M H+ , 
the supersymmetric loop form factors of these interactions are suppressed by one power of tan/3 
with respect to the terms described by C. On the contrary, this is no longer true for the Hb^bji 
and hbibn form factors [30]. For these couplings C is only correct in the limit Mg USY 3> M\. 



3 Quantum corrections to T(t -> bH + ), T(H + -+tb) 



The tree-level partial widths read 

•,2 



I*«(t-&17+) = ^|%| 2 m t 3 AV 2 (l, fe+ , ( 



(1 - Q H + + Qb) (cot 2 / 3 + % *an 2 /3) + Aq b ] , (22) 



(1 - r t - r b ) (r t cot 2 /3 + r b tan 2 /?) - 4 r t r b ] , (23) 



where we have defined the ratios — m 2 H+ /m 2 , — m b t /M 2 H+ and the A 1 / 2 term is a 

kinematic factor 

X(l,x,y) = l + x 2 + y 2 -2{x + y + xy). 

From now on, we shall assume \ V t b\ — 1 and neglect light fermion generations. For values of the 
parameter tan/3 ^ 15 (the inflexion point being given by tan/3 > ^Jm t /m b ~ 7) virtual quantum 
effects are largely dominated by the corrections to the right-handed bottom Yukawa coupling. In 
that limit the tree-level widths reduce to 

T^(t^bH + ) = _|!^|-(l-^ + ) 2 % tan 2 /?, (24) 

I*«(tf+-t6) = ^rM% + (l-r t ) 2 r b tan 2 /3, (25) 
in which we have also taken into account the smallness of m b as compared to m t , M H+ . 



3.1 Standard QCD correction 

As we have proved in Appendix B applying the OPE, both leading and subleading \og(Q/m b ) 
logarithms in the t — > bH + and H + tb renormalized decay widths can be resummed by using 
the running, a 2 corrected, bottom mass in the zeroth-order expressions. The one-loop finite QCD 
terms, though, are also sizeable, and have to be taken into account. In this section we derive 
improved expressions for the QCD-corrected decay rates, including both kind of effects, for large 
tan/3 values. 

The one-loop QCD-corrected expressions for the t (H + ) decay rates [18-21] can be greatly 
simplified after expanding them in a series in powers of r b (q b ) and retaining only the first-order 
term. As we are mainly interested in the region of large tan/3, we will provide formulae valid for 
those values of tan/3, for which eqs. (24), (25) apply. An explicit evaluation of the departure from 
this approximation for the one-loop result will be done in section 4. 
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In the H + case we perform a simultaneous expansion in powers of r b and r t . Retaining terms 
up to rf and considering the logarithmic factors to be of 0(1), the resulting approximation to the 
one-loop formula is 



1 + — 

7T 



° 2x M?,.i\ nrnunSix 



32vrM^ 



3 + 6 r t + r 2 t - ^ r\ + 2 log(r 6 ) + (-4 r t - y r\ - rf ) log(r 



(26) 



As can be seen from the above equation, there is no need to do the resummation of the \ogr t 
logarithms, as they are either small when r t is close to 1 or suppressed by at least a power of r t 
when it is small. 

In the limit of very small r t , eq. (26) reduces to 



a 



l + ^(3 + 21ogr 6 ) 



7T 



r (0) (#+ ->t&), 



(27) 



where we have introduced the quantity r(°), which is formally identical to T tree but has as input 
parameters the on-shell renormalized ones. The finite part in eq. (27), 3a s /ir, stands for a correction 
of about +10% (for a s ~ 0.1), whereas the full correction is large and negative, due to the much 
bigger logarithmic term. 

For the t — > b H + decay, the expansion in q b reads 



( t ^bH + ) ~ g2m * 



-^app 



(l-te+) 2 %tan 2 /? |l + ^ x 
2-5q H 



9 2tH 3 , q H+ , 

4" — + 2 log% -r^--: log ^ + + - % 



log(l-<z H+ ) 



+ log^+ log(l - +2ReLi 2 (l - g H+ ) 



(28) 



In the limit — ► 1, the ratio Tqcd/F^ becomes infinite and perturbation theory breaks 
down, as the 6-quark moves too slowly in the top rest frame. Nevertheless, the correction goes to 
zero due to the presence of the kinematic suppression factor. 

At this point we are ready to incorporate the resummation of the leading and next-to-leading 
q b , r b logarithms, as explained in section 2.2, which amounts to replacing m b in eqs. (26) and (28) 
by the running bottom mass at the proper scale. 4 The one-loop QCD-corrected widths are then, in 
the large tan/3 limit and including renormalization group effects up to next-to-leading order, given 
by the following improved (imp) formulae 



647TM2, 



1 + 



a s (m 2 ) 



7T 



7- — -2 log(l -q H+ )+2(l- q H+ ) 
+ Q + flog(l - q H+ )) (1 - q H+ f | , 



(29) 



4 We refer the reader to Appendix B for a proof of that statement. 
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rg&,(tf+-t5) = ^^M H+ {l-r t fml{Ml + )t^^ 

a 8 (M%+) ri7 , 16 , / 10 , 40 ., , , 

i + H -j +6rt+r * "27 1 + r * ~ y r * y n ' g(r * } 



(30) 



I 7r 

where a s (Q 2 ) is the MS-scheme running coupling constant and fn b (Q 2 ) the MS running mass 
expressed in terms of the bottom pole mass. 

Finite parts in Y mp and T app differ (see e.g. the 17/3 in eq. (30) and the 3 in eq. (26)). There 
is an implicit scheme conversion in going from eqs. (26), (28) to eqs. (29), (30): the bottom pole 
mass has been replaced for the running MS mass in the prefactor and the log(r fe ) has been absorbed 
into m b . Notice that the non- logarithmic terms of Y'qqd have been explicitly included in ^ 1 qcd-> as 
they are not accounted for by the renormalization group resummation techniques. 

3.2 Supersymmetric corrections 

The effective lagrangian prediction for the SUSY-QCD and SUSY-EW corrected decay rates can 
be read from eq. (21). No tan/3-enhanced vertex corrections contribute to the matching and the 
result is obtained by simply inserting the effective coupling, eq. (3), into the zeroth-order width 

?MsY = 1 ^2 r(0) • (31) 

SUSY (1 + AmJ 2 

We want to compare eq. (31) with the diagrammatic on-shell expressions for the one-loop SUSY- 
QCD and SUSY-EW corrected t^bH+,H+ -> t b partial widths [14, 16, 17], which we will denote 
by ^Isvsy- F° r l ar g e tan/3 values, the only sizeable diagrams are those that contribute to the scalar 
part of the bottom quark self-energy, entering the computation through the mass counterterm. For 
the SUSY-QCD corrections, the diagram that matters is shown in fig. 4. By simple power counting 
one can realize that it is finite. Moreover, neglecting O (m b / Mgy SY ) contributions, its value is 
essentially given by that of the three-point diagram in eq. (6): Am^ CT . 

Similarly, the diagram relevant to the SUSY-EW corrections is a two-point one with a chargino 
(neutralino) and a stop (sbottom) inside the loop. As for the SUSY-QCD case, it is finite, and 
its value can be approximated by the corresponding three-point diagram where an extra H® leg is 
attached to the scalar line. Its contribution is thus given by Am b EW in eq. (9). 

Collecting the results from eqs. (6) and (9) via eq. (5), the one-loop SUSY corrected decay rates 
can be cast into the formula 

r'susr = (l - 2 Am,) r<°> + AT susy ■ (32) 

The term AT susy, which contains non-universal and tan/3-suppressed contributions to the decay, 
is very small provided tan/3 is large, as we have numerically checked. 

Both prescriptions, eqs. (32) and (31), are equivalent at first order in perturbation theory (PT) 
and consequently do not differ significantly when the corrections are small. In general, though, 
Am b can be a quantity of 0(1) for large enough tan/3 values, in which case eq. (31) is preferred as 
it correctly encodes all higher-order Am b effects (see the discussion in section 2.1 and in Appendix 
B). 

3.3 Full MSSM renormalization group improved correction 

In section 2.2, we saw how the effective lagrangian (21) accounts for the higher-order tan/3-enhanced 
SUSY quantum corrections and also for the leading and next-to-leading QCD logarithms, including 
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those in diagrams like fig. 2. We define the improved values for the decay rates of the two processes 
under study in the MSSM as 

1 



-pimp 



-,imp 



MSSM — L QCD 



+ AT SUSY • 



(1 + Am b ) 

which also incorporates the one-loop finite QCD effects. Neglecting the small tan/3-suppressed 
ATsusy effect, one has 



9 2 N C 



H+< 



(l + Am b Y 



tan 2 /? x 



« s (m 2 + ; 



7T 



17 2 16 o / A 10 2 40 , , , 

— +§r t +r t - —r t + l-Ar t - —r t - —r t ) log(r t 



(33) 



■pimp 
1 MSSM\ 



■2 m 2 b (m 2 t ) 



(t^bH + ) = — ~r~T7T~ m t (1 — Qh+) 9 tan 7? 

V ' 64vrM^ 1 y qH+ > (1 + Am b ) 2 



1 , 



7T 



7 " IT " 2 log(1 " ^ +) + 2 (1 " ^ +) 



+ Q + ^iog(i-^ + ))(i-^+) 2 } • 



(34) 



The above formulae contained all the improvements discussed in this article. In order to compare 
them to the diagrammatic one-loop MSSM results, we introduce r^fgsjj^ 



pi— loop 

^1—loop _ -pimp 1 SUSY 
MSSM — 1 QCD -p(o) ' 



(35) 



which only differs from ^mssm m ^ na t no resummation of the SUSY-QCD, SUSY-EW corrections is 
performed. Comparing ^ms'sm-i ^mssm one can assess the size of the higher-order tan/3-enhanced 
effects. 



4 Results on the decay rates 

Although the t — > bH + and H + — ► tb decays are mutually exclusive, in the effective 2HDM 
lagrangian we constructed in section 2, the supersymmetric corrections to both observables are 
encoded in the same effective coupling. Therefore, we prefer to present the study of these corrections 
simultaneously, stressing the points they have in common. 

To quantify the importance of the quantum corrections we introduce the relative correction to 
the width S, defined as 

p -ptree 

ST X = -^—^ . (36) 

4.1 One-loop vs. NLO-improved QCD corrections 

Figures 5 and 6 analyse the gluonic corrections to the t — > b H + and H + — ► tb decay rate, showing 
their dependence on the mass of the charged Higgs boson for tan/3 = 10 and 30. The dotted lines 
represent the relative shifts, STqq°^ p (36), produced by the one-loop QCD corrections, ^qc°£) P , 
which have been computed using the formulae in refs. [19, 21]. 
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In the limit of large tan/3 and small q, , r b , the above one-loop results admit simpler approximate 
expressions, which we have derived in eqs. (26), (28). These approximations have a lower bound of 
validity, which can be roughly set at tan/3 = 10. In this paper we will not consider values of tan/3 
smaller than 10. Inserting the expansions (26), (28) into eq. (36), one obtains the corresponding 
relative shifts SVq^ d for the t — > bH + and H + — > tb partial widths, which are plotted using the 
dashed lines in figs. 5 and 6 respectively. 

For tan/3 = 10, the first term in the r b expansion of the non tan/3-enhanced one-loop QCD 
corrections to the t — ► bH + decay rate stands for a contribution of about 5%. For the sake of 
simplicity, we omitted this term in eq. (28), but we have included it when drawing the ST ( qq D 
curve in fig. 5. The extra correction is almost negligible for the tan/3 = 30 curve. In the H + — ► tb 
decay rate, fig. 6, and for tan/3 > 10, eq. (26) is always extremely close to the one-loop result, and 
the ST'qq D curves are not shown. 

As can be seen in fig. 5, a discrepancy appears between ST^q D and STqqq P close to the 
threshold, which can be traced back to the fact that we dropped the m b kinetic terms in the 
approximated formula. Similar problems should be present in the H + — > tb case, fig. 6, when 
approaching the threshold, but our plot starts at a conservative M H+ = 250 GeV value for which 
the truncated series, eq. (26), with r t = 0.5, is still valid. In any case it makes no sense to try to 
include higher-order r™ terms because close to the threshold the perturbative expansion is no longer 
reliable: the decay products move slowly in the decay particle's rest frame, and long-distance non- 
perturbative effects can significantly modify the perturbative prediction. Moreover, in this region 
the branching ratio is very small and therefore the corresponding decay channel loses its relevance 
for the charged Higgs phenomenology. 

As was justified in section 2 using the operator product expansion, the replacement of the renor- 
malized bottom mass and strong coupling by their running two-loop MS values correctly resums 
leading and next-to-leading r b , q b logarithms. In eqs. (29), (30) the substitution was explicitly 
done. In fig. 5 the numerical effect of the improvement corresponds to the difference between the 
dashed (ST'qq D ) and solid (ST 1 ^q D ) curves. For the H + — ► tb decay the improvement is essentially 

given by the difference between the dotted (^qc°d P ) an d solid {S^qqjj) curves. 

Even for moderate tan/3 values around 10, the QCD corrections are larger than 50%, driven 
by the big q b , r b logarithms. The resummation of the leading logarithms is mandatory, specially 
for the H + decay where logr^ is unbounded as M H+ increases. The effect of the LO and NLO 
resummation diminishes the top partial decay rate in about 5% and the charged Higgs decay rate 
in about 15%. 



4.2 Supersymmetric corrections 

Figure 7 focuses on the genuine supersymmetric corrections to the H + — > tb partial width. As they 
are dominated by the universal Am b effect, the results for the H + — > tb plot represent fairly well 
the effects of the corrections on t — ► b H + too. Curves are shown for two values of the /x-parameter 
and for two different sparticle spectra. 

In the "heavy" spectrum, the gluino and the lightest sbottom and stop have a common 1 TeV 
mass. The squarks and gluinos are nearly degenerate and they are much heavier than the H + mass, 
justifying the use of the effective lagrangian approach. As only tan/3 > 10 values are considered, the 
approximation consisting in neglecting the non-universal and tan/3-suppressed terms denoted by 
AT susy in eq. (32), which is represented by the dashed ^susy curves, fits very well the one-loop 
calculation (the latter is not shown in this case). The corresponding effective lagrangian prediction, 
eq. (31), which includes all tan/3-leading terms appearing at higher orders in PT, is represented by 
the solid 5T e M QV lines. 
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t-> H + b 
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m H+ (GeV) 

Figure 5: Comparison of the QCD contributions to the t — > 6^T + decay width, as a function of 
M H+ , for tan/3 = 10 (upper set) and 30 (lower set). The dotted line denotes the one- loop result [21], 
the dashed line the approximation of eq. (28), and the solid line the NLO-improved one in eq. (29). 
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Figure 6: Comparison of the QCD contributions to the H + tb decay width, as a function of 
M H+ , for tan/3 = 10 (upper set) and 30 (lower set). The dotted line corresponds to the one-loop [19] 
correction, and the solid line to the NLO-improved result, eq. (30). 
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A second, lighter, spectrum is denned by M~ g = 500 GeV, and the masses of the lightest sbottom 

and stop around 200 GeV. The curve labelled ST^}g Y corresponds to the full one-loop computation, 
including all possible gluino, chargino and neutralino loops. Even for this light spectrum and for the 
chosen set of parameters, 5T a g^ SY gives a good estimate of the one-loop correction. This illustrates 
the fact that our effective lagrangian C in eq. (21) describes the charged Higgs interaction correctly 
even if M SUSY < M R+ . It shows that Am b accounts for most of the effects and we can trust the 
validity of the improved result. 

Typical values we found for the SUSY correction are 15%-30% with the heavy spectrum and 
~ 40% with the light one. In both cases, the results depend heavily on the ji and tan/3 parameters, 
the size of the correction growing almost linearly with their absolute values. Although not shown 
in the plots, the main contribution to ST SUSY comes from the SUSY-QCD diagrams. Only for a 
very large At values can the electroweak corrections be comparable. 

The SF'susy curves correspond to the relative correction to the widths as evaluated using 
eq. (31), an expression derived from the effective lagrangian in section 2. While <5T^°?, 5T a <$ SY 
do not include higher-order Amj) effects (which can be potentially of 0(1)) these tan/3-dominant 
effects are correctly resummed to all orders in PT in the expression for 8Ygjj SY . 

The difference between ST e JjJ SY and S^'susy ^ rs ^ appears at order (AmJ 2 , and is always 
positive, opposite to the negative standard QCD corrections, for Am b > —1.5. 5 Therefore, for 
negative (positive) values of Am b , that is, positive (negative) corrections ^r^^y, the higher-order 
terms tend to reinforce (suppress) the correction. As Am b is mainly given by the SUSY-QCD 
contribution, eq. (6), this correlation is seen in association with the sign of \x. 

Just to give some examples, for a negative S^ l sv°SY = 30% correction, which corresponds to 
Am b = —0.15, the extra higher-order terms contained in ST e s {j SY increase the partial width by 8%. 
For Am b = —0.2, a number that can be obtained from eq. (8) by setting tan/3 = 20, a s = 0.1, the 
difference between ST e s ^j SY and S^sjjsy ^ s °f or der +16%. 

The only restriction to the potential size of 5T SUSY is set by the renormalized bottom Yukawa 
coupling, which is required to remain perturbative from the GUT scale to the scale of the corre- 
sponding decay. This is guaranteed in our calculations by demanding < 1.2 at low energies (see 
e.g. [13]), implying the following combined bound on tan/3 and Am b : 

Am b >-^r 9 m ^ m ^ tan/3 - 1 ~ 0.014 tan/3 - 1 . (37) 
1.2 y/2M w 

In the above example, with tan/? = 20, the minimum allowed value for Am b is —0.72. If eq. (8) 
holds for negative fi, and using ot s {M SUSY ) ~ 0.1, it is found that a maximum allowed correction, 
6T< SUSY ~ +200%, is obtained around tan/3 = 40. 



4.3 Full MSSM correction 

We shall now show the combined effects of the QCD, SUSY-QCD and SUSY-EW corrections in 
the partial decay widths under study, starting from three different sets of curves: ST l ^ D , i.e. the 
QCD correction including the renormalization group resummation of the bottom mass logarithms 
up to NLO; ^ms'sm^ the full one-loop MSSM contribution as defined in eq. (35); and the MSSM- 
improved contribution, ST l ^g SM , defined in eq. (34). 

Figure 8 shows the dependence of the relative corrections to the width 5F(t — > bH + ) on the 
mass scale M SUSY , defined as a common value for the gaugino mass, M2, the gluino mass and the 

5 The comparison is between 1/(1 + Am b ) 2 and 1 — 2Am b , eq. (32), the approximated one-loop result as denned 
in this paper and in [17, 15]. 
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H + -> t b 




Figure 7: The SUSY contributions to the H + — > tb partial decay width, as a function of tan/3, for 
M H+ = 350 GeV and two values of fi. The dashed lines denote the approximation ATsusy = of 
eq. (32), whereas the solid lines correspond to the effective width, eq. (31). For the heavy spectrum 
one has M~ g = = = 1 TeV, b±, t\ being the lightest sbottom and stop respectively. A t = 
500 GeV, the \i values are shown in the plot. For the light spectrum we have set M~ = 500 GeV, 
mi = 250 GeV and m^ = 180 GeV. In this case, we also show a dotted curve corresponding to 
the one-loop result of ref. [15]. 
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masses of the lightest stop and sbottom. As we keep the value of \x fixed, the SUSY contribution 
smoothly goes to zero like ^/M SUSY when M SUSY increases. Contrarily, if all mass parameters 
are sent to infinity together, the SUSY correction tends towards a constant value, determined by 
Am b ~ ±(a s /37r) tan/3, eq. (8). A similar behaviour occurs for 5T(H + — > tb) with a different 
renormalized value for 5Tq CD . 

The difference between 5T l q^ D and S^mssm is due to the SUSY corrections, which were already 

considered in the above section. The mismatch between ^^ms'sm ano - ^ % mssm ls produced by the 
tan/3-enhanced higher-order effects that are resummed in the latter. 

Figure 9 shows how the full MSSM correction evolves with tan/3. While STq CD has a mild 
dependence on tan/3 that is almost saturated around tan/3 = 20, the SUSY part gets more and 
more important as tan/3 increases. One can see that for the chosen parameters <5r 



SUSY 



becomes of 

0(10%) around tan/3=30. For negative values of fi, of 0(M SUSY ), and for sufficiently large tan/3 
values, the total correction can be considerably reduced with respect to the naive QCD prediction. 
A similar behaviour is found for 5T(H + — ► tb). 
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Figure 8: Evolution of the corrections to the t — ► bH + width, for M H+ = 125 GeV and tan/? = 30, 
as a function of a "common SUSY mass" , M SUSY = M2 = M~ = = , and A t = 500 GeV. 
The dashed line corresponds to the QCD-improved width, eq. (29), the dotted line denotes the 
one- loop MSSM result, eq. (35), and the solid line denotes the MSSM-improved one, eq. (34). 
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Figure 9: The corrections to the t — > bH + width for M H+ = 125 GeV as a function of tan/3. The 
rest of the parameters are those of the heavy spectrum in fig. 7. The dashed line corresponds to 
the QCD-improved width, eq. (29), the dotted line denotes the one- loop MSSM result, eq. (35), 
and the solid line denotes the MSSM- improved one, eq. (34). 
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5 Results on the branching ratios 



Above we have described the effects of the QCD, SUSY-QCD and SUSY-EW corrections on the 
decay widths of t — ► b H + and H + — > t b as a function of the MSSM parameter space. In the case 
of t — ► bH + , assuming that the only other possible decay channel is t — > bW + , we shall present 
the results on the BlZ(t — > bH + ) and we shall use these computations to exemplify how much the 
radiative corrections implemented here can change the actual reach of the Tevatron collider in the 
search of H + in the indirect mode, missing leptons/dileptons in the t — > 61U + decay. 

Here, the results from the frequentist analysis of DO indirect H + searches [5] are used to derive 
constraints on the tan/3-M^ + plane (see e.g. [6] for results on similar indirect H + searches by the 
CDF collaboration). 

In fig. 10 we draw curves of constant BlZ{t — > bH + ) based on ^QQ D {t — > bH + ), eq. (29), and 
including the one-loop QCD corrections into the computation of T(t — > 6IU + ). We do not show 
curves that have a branching ratio smaller than 0.1 because, for such regions of parameters, the 
t — > &.£f + decay channel has little phenomenological relevance. The grey area at the bottom-right 
corner of the figure is the region excluded by the DO frequentist analysis data. 

The plots in fig. 10 compare to the plots in fig. 11. Here we show curves of constant BTZ(t — ► 
bH + ), using the MSSM-improved formulae for the partial t — ► bH + decay rate, eq. (34). The soft 
SUSY-breaking masses are chosen to produce a heavy SUSY spectrum, with M~ g = = m~ = 
1 TeV. As in fig. 10, the dark area on the bottom-right corner corresponds to the experimentally 
excluded region. 

For positive values of Am b (left plot in fig. 11), both QCD and SUSY-QCD corrections reduce 
the tree- level partial width of the t — > b H + decay channel, and the bound on the BIZ moves to higher 
tan/3 values. In our example plot, with fi = 500 GeV, the excluded region starts at tan/3 > 100 
and it is not shown. Conversely, for negative Am b values, the supersymmetric corrections partly 
compensate for the QCD reduction of the width, and the bound is found for lower tan/3 values. 
This fact can be checked in the plot on the right of fig. 11, corresponding to fx = —500 GeV. Values 
larger than 0.4 for BlZ(t — > bH + ) are obtained when tan/3 <; 55. The experimental bound starts 
around tan/3 = 65, in a region where hb(m t ) > 1.2, which implies that the bottom Yukawa coupling 
becomes non-perturbative below the GUT scale [13]. This fact is denoted in the plots by changing 
from solid to dashed line style. The same remark applies for fig. 10. 

The H + — > tb branching ratio, which is expected to be tested at the LHC and at the NLC, 
is depicted in fig. 12. On the left plot, contour lines of constant BIZ are drawn using the QCD 
improved width, eq. (30). Similarly, the right plot shows curves of constant B7Z{H + — ► tb) for the 
MSSM-improved result, eq. (33), with /j, = 500 GeV, the rest of SUSY parameters being equal to 
those of fig. 11. It has been assumed that no decays of H + into pairs of R-odd SUSY particles [7] 
were possible. This is guaranteed by the choice of the soft SUSY-breaking masses and by cutting 
the plots at M R+ = 500 GeV. 
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Figure 10: Curves of constant branching ratio for the t — > bH + channel. The figure shows the 
QCD-improved, eq. (29), result. The transition between the solid and dashed styles occurs when 
the bottom Yukawa coupling crosses the bound hb(m t ) < 1.2. As explained in the text, this bound 
guarantees the perturbativity of the Yukawa up to the GUT scale. Finally, the shaded area defines 
the 95% C.L. exclusion boundary in the t&nf3-M H+ plane for m t = 175 GeV and a(tt) = 5.5 pb 
that can be derived from the DO frequentist analysis in ref. [5]. 




Figure 11: As fig. 10, but plotting the MSSM-improved result, eq. (34), for fi = 500 GeV (left plot) 
and [i = —500 GeV (right plot). The rest of relevant SUSY parameters are given by M~ g = M2 = 
= m~ b = 1 TeV, A t = 500 GeV. In the fi = —500 GeV plot, the shaded area is excluded by the 
DO frequentist analysis in ref. [5]. 
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Figure 12: Curves of constant branching ratio for the H + tb channel. On the left, the QCD- 
improved values, eq. (30), on the right, the MSSM-improved result, eq. (33). The parameters 
chosen for these plots are fj, = 500 GeV, M~ = M 2 = m f = m~ = 1 TeV, A t = -500 GeV. 
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6 Conclusions 



Using an effective lagrangian description of the MSSM, we have investigated the virtual super- 
symmetric effects that modify the tree-level relation between the bottom Yukawa coupling and the 
bottom mass, which are dominant in the large tan/3 regime. Motivated by the fact that these effects 
do not vanish for large values of the SUSY masses and are potentially of 0(1), we have derived the 
expressions for the bottom Yukawa couplings that resum all higher-order tan/3-enhanced quantum 
effects. These expressions have a natural interpretation and are easily deduced in the context of 
the effective lagrangian formulation. We have also shown that they can be equivalently deduced in 
the framework of the full MSSM. 

As an interesting application of our results, we have computed the partial decay rates for the 
t — > bH + and H + — > tb decay channels, relevant to supersymmetric charged Higgs searches at 
present and future colliders. First we have considered the QCD quantum corrections to these 
processes and, applying the OPE, we have performed the resummation of the leading and next- 
to-leading logarithms of the form log Q/m b . Concerning the supersymmetric corrections, we have 
compared our results with those of previous diagrammatic one-loop analyses in the literature and 
we have shown the numerical relevance of the resummation of the tan/3-enhanced effects derived 
in this work. Collecting the above improvements, we have finally computed the corresponding 
branching ratios, BlZ(t — ► b H + ) and B1Z(H + — > tb). As an example, we have shown, for different 
sets of the MSSM parameters, the effect of the quantum corrections in determining the region of 
the M H+ -tanl3 plane excluded by the DO indirect searches for a supersymmetric charged Higgs 
boson in the decay of the top quark. 
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Appendix A The effect of Am b at all orders 



In this appendix we perform the resummation, in perturbation theory, of the leading supersym- 
metric effects contained in Am b and find agreement with the effective lagrangian result of sect. 2. 

When relating tan/? to physical decay rates or to quark masses and Yukawa couplings one 
also has to address the question of the proper treatment of standard QCD corrections. This 
issue seems to come into play in the very beginning when one defines the renormalization scheme 
and scale for the bottom quark mass, which, for example, enters the off-diagonal elements of 
the 6-squark mass matrix. Here we want to stress that the issues of large supersymmetric tan/3- 
enhanced corrections related to the diagram in fig. 4 and the treatment of standard QCD corrections 
related to gluonic corrections to the quark self-energy can be treated independently of each other. 
In the following we shall concentrate on the tan/3-enhanced SUSY-QCD corrections, induced by 
supersymmetric particle loop effects. After performing the resummation of these corrections to all 
orders in perturbation theory, one can subsequently include the standard QCD corrections, whose 
proper treatment is discussed in sect. 2. 

The quantity Am b is proportional to (a s / 'it) (fi tan /3/M SUSY ) and of 0(1) when, simultaneously, 
/i ~ M SUSY and tan/3 is large. In that case one should resum its effects to all orders in PT to 
obtain a reliable prediction. As was shown in section 2, the first thing one should realize is that 
there are no higher-loop diagrams contributing to the mass renormalization (nor to the decay rate) 
of order a"tan n /3 with n > 1. Diagrams with extra /xtan/3 insertions are suppressed by powers of 
m b/MsusY- This can be easily seen in the effective lagrangian approach, where such contributions 
would arise from higher-dimensional operators with more Higgs boson fields, whose couplings are 
suppressed by extra powers of M SUSY - 

Different renormalization schemes use different values for the renormalized bottom Yukawa 
coupling hb [22]. In theories with spontaneous symmetry breaking, though, there is always a link 
between the value of hb and the physical bottom mass, m b : the dressed bottom propagator must 
have a pole for on-shell external legs, or conversely the inverse propagator must vanish. At one 
loop this relation reads, considering only the gluino corrections 

hbVi + 5hbV i + hbVi Am b = m b , (38) 

5hb being the counterterm of hb- The l.h.s. of the previous equation is graphically depicted in 
fig. 13. 

9 




Figure 13: Feynman diagrams contributing to the bottom pole mass up to first order in PT. From 
left to right, the renormalized bottom mass, the bottom mass counterterm and the finite one-loop 
Feynman graph contributions are shown. The dashed line in the last diagram denotes a sbottom 
and the solid line a gluino. The cross represents the insertion of the bottom mass counterterm. 

We are not displaying the wave function renormalization to avoid an unnecessary complication 
of the argument. Note that v\ receives no one-loop QCD corrections and thus its renormalization 
only adds effects suppressed by aE\v/oi s , which allows us to identify 8hbV± with Sm b . Besides, in 
any renormalization scheme one has = hbV\, with and hb denoting renormalized quantities. 
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Therefore, 

(h b + Sh b ) v 1 = m b H + 5m b , 



(39) 



and one obtains, at first order 

m b + 5m b = m b — m b Am b . (40) 

The l.h.s. of eq. (40) is just the bare bottom mass, m b . 

When evaluated beyond first order, scheme differences appear in the r.h.s. of eq. (40). In the 
on-shell scheme, the renormalization condition being given by m b = h b v\ = m^, one would obtain 
that the bare bottom mass is equal to m b (1 — AmJ, while in the MS-scheme, for which 5m b is 
zero as Am b is finite, one would have m b /(l + Am b ). Both results are equivalent at first order in 
Am b , as they should. 

To proceed with the resummation, we come back to the relation between the Yukawa and the 
pole mass. Although no n-loop diagrams produce a™ tan n /3 corrections for n > 2, there is one and 
only one genuine n-th-order diagram left (see fig. 14), which contains the insertion of a (n — l)-loop 
counterterm into a one-loop diagram. Then, all dominant terms in the large tan/3 limit, at all 
orders in PT, are contained in the equation 

h b v\ + 5h b v\ + h b v\ Am b + 8h b v\ Am b = m b . (41) 




Figure 14: Full set of SUSY-QCD dominant diagrams, in the large tan/3 limit, contributing to the 
bottom pole mass at all orders in PT. The first three diagrams are those of fig. 13. In the fourth 
one, the cross denotes the insertion of the Sh b counterterm, and the solid and dashed lines denote 
gluino and sbottom propagators respectively (see fig. 16). 

Beyond tree level, the bLbnH® coupling is no longer equal to h b , so it is denoted by h b , with 
counterterm Sh b . 6 This fact was not important in eq. (38) because we were just considering the 
first-order result. 

Before proceeding, one technical point in equation (41) deserves further clarification. The last 
term in the l.h.s. corresponds to the true three-point diagram in fig. 14. In the large tan/3 limit, 
though, its value, Sh b viAm b , coincides with the two-point contribution, h b viAm b , after replacing 
the renormalized coupling by the counterterm. A derivation of this result is written at the end of 
this appendix. 

The last step in our argument is to justify the equality 

h b + 5h b = h b + 5h b , (42) 

which can be regarded as the identity of the bare quark and squark Yukawa couplings, which is 
guaranteed by the underlying supersymmetry governing the relations between the bare lagrangian 
parameters in the ultraviolet. 7 No soft SUSY-breaking dimensionful couplings can induce modifi- 
cations to eq. (42), allowing for the extraction of a common h b v\ + 5h b v\ factor in (41). At the level 

6 The tree-level coupling is in fact m b fj, tan/3, but again neither fi nor tan/3 receive QCD corrections at first order. 
7 This is true if a regularization method preserving SUSY is used, such as dimensional reduction. Deviations from 
eq. (42) in the MS-scheme will be loop-suppressed, not affecting the conclusions of this appendix. 
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of bare couplings one does not need to make reference to any particular renormalization scheme. 
Thus, one has 



(h b + Sh b ) v 1 = m b H + 6m b = b , (43) 

1 + Afiij 

where the r.h.s. is expressed in terms of physical quantities, Am b being independent of m b . 

For the rest of this appendix we will derive expressions valid to all orders in PT in the large tan/3 
limit for the H + and H, h, A dressed couplings to tb and bb respectively, recovering the effective 
lagrangian results one can find in [24]. The calculation involves contributions from three-point loop 
diagrams with one external on-shell Higgs leg whose momentum we have neglected. In section 4, 
the departure from this assumption for the H + and t decay rates has been shown to be small, as the 
extra contribution inducing the momentum dependence does not include any tan/3 enhancement 
factor. More complete formulae including the momentum dependence for the decay rates of the 
neutral Higgs bosons can be found in ref. [30]. 

Let us start with the simplest case, that of the charged Higgs H + and of the pseudoscalar A, 
for which there are no vertex loop diagrams tan/3-enhanced with respect to the tree- level coupling. 
The relevant Feynman diagrams are just the tree-level Yukawa and the counterterm. From eq. (43), 
the renormalized decay amplitudes are given by 



i (h b + 6h b ) sin PH + iP R b = , tan/3 H+tp b 

1 + Am, v 





., rj . sin/3 - rattan/? - 

. iht + SI%) _ A , mb _ __^_ J 4 675i> . (44) 

Therefore, in this case, the result of the resummation is to effectively modify the tree-level Yukawa 
coupling by the universal 1/(1 + AmJ factor. 

The case of the CP-even neutral Higgs bosons is a little bit more involved. Depending on the 
relation between (5 and the mixing angle, a, the one-loop correction to the vertex diagrams can be 
importantly enhanced. The full set of potentially relevant graphs is shown in fig. 15. 

k 7 bb 

+ < kg + K 





b a V b a 



Figure 15: Vertex diagrams contributing to the renormalization of the Higgs-fermion Yukawa inter- 
action. From left to right, the renormalized Yukawa coupling, the Yukawa counterterm, the one-loop 
contribution and the higher-order diagram containing the insertion of the 5h b counterterm. The 
solid and dashed lines inside the loops denote gluino and sbottom propagators respectively. The 
cross in the fourth diagram denotes the Sh b counterterm. 

One obtains, for the Hbb renormalized amplitude 

.,, . cos a -, ~ \ sinaAm ft ttT7 

-i (h b + 5h b ) — pr- Hbb -i[h b + 5h h ) —= £ Hbb = 

a /9 V / a tan H 



s/2 v > V2 tan/3 

m h cos a f tan a \ 

r- , b — 1 + Am b 

V2(l + Am b )vi V tan/3 / 



rriicosa / tana\ - , N 

-i b — 1 + Am b Hbb . (45) 
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Again, the resummation amounts to the inclusion of the universal 1/(1 + AmJ factor. However, 
there is an additional Am b term inside the parenthesis, which constitutes the non tan/3-suppressed 
contribution coming from the SUSY-QCD vertex diagrams. Similarly, for the hbb one has 



. . sma - . /- ~ \ cos a L ^ m b , t, 

i (rib + ohb) — ■=- hbb — i [fib + 0/1& — ■= hbb 

V2 v ' V2 tan/3 

m^sina / Am, 



y/2 (1 + Am b ) vi V tan a tan/3 



hbb . (46) 



It can easily be checked that for large M H+ values, the limit that corresponds to the effective 
decoupling of one of the Higgs doublets, one recovers the SM hbb coupling 

whereas the Hbb coupling, H being heavy, still "feels" the decoupled sector 
m b tan/3 



y/2 (1 + Am b ) 



Hbb. 




Figure 16: The fourth self-energy diagram in fig. 14, shown in greater detail. A gluino propagator 
is denoted by the solid line inside the loop. The dashed lines denote sbottom propagators and the 
cross the insertion of a Shb counterterm. 



Two-point-three-point diagram identity 

Let us evaluate the amplitude associated to the three-point Feynman diagram of fig. 16. Neglecting 
the external momentum, it can be written 

-(8na s ) C F M^tau/3 J — (fe2 _ ^ _ _ bP L b , (47) 

where Cir = 4/3isa colour factor and the two-dimensional rotation matrices Z transform the weak 
eigenstate sbottom basis into the mass eigenstate basis. Expressed in terms of the mixing angle 9 b , 
the components of Z read: Z\\ = Z22 = cos^, Z\i = —Z21 = sin^^. The term between parentheses 
in the numerator of (47) and the combination ShbV2 = ShbVi tan/3 come from the counterterm to 
the tan/3-dominant interaction H^bfjpL, after the Higgs field develops its vacuum expectation value 

V2- 
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Splitting the implicit i,j sum into the i = j part and the rest of the terms we obtain 

k . 2 f d n k 1 / 1 1 \ - 

2 Sm b J (2^T fc 2 - Mf (^-m^) 2 + (fc 2 -m? 2 ) 2 J 

/d n k 1 
(^ (^-M 2 )^-^)^ 2 -^) ^' (48) 

the constant k being a short-hand for the constant prefactor of the integral in (47). 

The second term in (48) is of the same form as Am b . Adding and removing this term times 
tan 2 29 b and rearranging terms one arrives at 

f d n k 1 - 

K J (2tt)" (k 2 - M?)(k 2 - m? )(fc 2 - m 2 J L 

k 2 f d n k 1 ( m i ~ m i ) 2 

+ j sin22 ^y (2vr)« k 2 - M? (k 2 - m 2 bi ) 2 {k 2 - mff ^ (49) 

Now one can make use of the tree-level relation 

= H^-z'M (50) 

mf — mf 

01 02 

to write (A is dropped since it is not tan/3-enhanced) 
d n k 1 



(2tt)" (A; 2 - M~ 2 )(fc 2 - m- 2 )(fc 2 - m? 



/cpfc 1 
(2^(fc 2 -M|)(fe 2 -m 6 2 i ) 2 (fc 2 -m? 2 ) 2 ^ 6 - (51) 

The second integral in (51) has two extra propagators and thus in the limit of heavy SUSY 
masses it is of 0(l/Mg USY ), whereas the first one is of 0(l/Mg USY ). One can conclude that the 
two- and three-point loop diagrams in fig. 14 are just related by 5hb/h , apart from contributions 
that are suppressed by powers of either tan/3 or m 2 tan 2 f3/Mg USY . The amplitude for the diagram 
in fig. 16 reduces to 



i ShbVi Am b 



b 

M'iusY " ~ V M 2 SUSY J ' - V tan/3 



bP L b. (52) 
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Appendix B Large logarithms in decay rates 



Our effective lagrangian C in eq. (21) contains the large logarithms associated with the running of 
the Yukawa couplings to all orders in perturbation theory. In general this procedure does not sum 
all the large logarithms that appear in a specific cross section or decay rate. In this appendix we 
show that for T(t — > b H + ) and T(H + tb) such additional, process-specific logarithms do not 
occur except in highly power-suppressed, numerically negligible terms. 

Let us first consider the decay H + — > tb: the optical theorem relates the decay rate to the 
imaginary part of the H + self-energy: 



T(H + ^tb) = -J— Imt [ d 4 xe- iq - x (H + \Tj\x)J(0)\H^ 
M H+ J 



(53) 



Here 



J(x) = 9 -^ff V tb tan/3 H+t L b R (x,Q) 
V2M W 1 + Am b 

is the scalar current stemming from the Yukawa interaction in eq. (21). All currents and couplings in 
this appendix are considered to be renormalized using a mass-independent renormalization scheme 
such us the MS scheme [28] . For the moment we also assume this for the quark masses and discuss 
the use of the pole mass definition, which is commonly used for the top mass, later. The decay rate 
involves highly separated mass scales m b <C M H+ , m t . First we assume that M H+ and m t are of 
similar size so that log M H+ /m t is not dangerously large. We return to the case m b <C m t <C M H+ 
later. To prepare the resummation of the large logarithm log m b /M H+ , we first perform an operator 
product expansion of the bilocal forward scattering operator in eq. (53): 

i fd 4 xe- i i- x (H + \Tj\x)J(0)\H + ) = C n (q 2 , m t , q) ( H + \O n \ H + ) (m b , Q) . (54) 

Here all dependence on the heavy mass scales m t and q 2 = Mjj + is contained in the Wilson 
coefficient C n , while the dependence on the light scale m b resides in the matrix element of the local 
operator O n . Both depend on the renormalization scale Q at which the OPE is carried out (so 
that Q is sometimes called factorization scale). The OPE provides an expansion of T{H + — ► tb) 
in terms of (m b /M H+ ) 2 . Increasing powers of m b /M H+ correspond to increasing twists of the local 
operator O n . Here the twist is defined as the dimension of the operator O n minus the number of 
derivatives acting on the Higgs fields in O n . 

The OPE in eq. (54) is depicted in fig. 17 where also the leading twist operator 0\ = m 2 (Q)H + H~ 
is shown. At leading twist the OPE, depicted in fig. 17, is trivial: the matrix element ( H + \0± \ H + ) 
simply equals m 2 (Q) and the Wilson coefficient C\ can be read off from eq. (33). In the leading 
order (LO) of QCD it reads 

ImCi = ^hrM 2 H+ (l~r t ) 2 1 tan 2 /3. (55) 

The QCD radiative corrections in T contain powers of the large logarithm a s log m b /M H+ . The 
OPE in eq. (54) splits this logarithm into a s logQ/M^ + + a s \ogm b /Q: the former term resides in 
the coefficient function C\ while the latter is contained in the matrix element ( H + \0± \ H + ). If we 
choose Q = 0(m t , M H+ ), then the logarithms in the Wilson coefficient are small and perturbative, 
but \ogm b /Q in the matrix element is big and needs to be resummed to all orders. One could 
likewise choose Q ~ m b and resum the large logarithm in the Wilson coefficient, but the former 
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Figure 17: The OPE in (54) to leading order in m b /M H+ and a s . The self-energy diagram on the 
left represents the left-hand side of eq. (54). The right diagram depicts (H + \0\ \ H + ). 

way is much easier here. In order to sum \ogm b /Q we have to solve the renormalization group 
(RG) equation for 0\. Since the Higgs fields in 0\ have no QCD interaction, the solution of the 
RG equation simply amounts to the use of the well-known result for the running quark mass m b (Q) 
(see eq. (18)) at the scale Q = 0(m t ,M H+ ) in 0\. In the next-to-leading (NLO) order one has 
to include the 0(a s ) corrections to T in eq. (33). First there are no explicit one-loop corrections 
to ( H + 1 0\ | H + ) , so that in the NLO we obtain \m.C\(Q) by simply multiplying the result in 
eq. (55) with the curly bracket in (33). Secondly in the NLO we have to use the two- loop formula 
for m,b{Q) in the matrix element. Since one is equally entitled to use Q = M R+ (as chosen in (33)) 
or Q = m t or any other scale of order m t , M H+ , there is a residual scale uncertainty. This feature 
is familiar from all other RG-improved observables. To the calculated order a s this uncertainty 
cancels, because there is an explicit term a s \ogQ/M H+ in the one-loop correction, so that the scale 
uncertainty is always of the order of the next uncalculated term. In our case this is 0(a 2 s ) and 
numerically tiny. In conclusion, our OPE analysis shows that at leading order in m b /M H+ all large 
logarithms in T(H + tb) can indeed be absorbed into the running quark mass in our effective 
lagrangian in eq. (21). Some clarifying points are in order: 

1) The summation of large logarithms in the NLO does not require the calculation of the two- 
loop diagrams obtained by dressing the diagram in fig. 17 with an extra gluon, as performed 
in [31]. This calculation only gives redundant information, already contained in the known 
two-loop formula for the running quark mass. 

2) At the next-to-leading order the result depends on the chosen renormalization scheme. Chang- 
ing the scheme modifies the constant term 17/3 in eq. (33). After inserting the NLO (two-loop) 
solution (18) for the running mass, this scheme dependence cancels between this term and 
j(f) 

in eq. (19). In the literature, sometimes, the one-loop result for T is incorrectly combined 
with the one-loop running bottom mass resulting in a scheme-dependent expression. 
No running top-quark mass is needed for the case m t ~ M H+ , and one can adopt the pole 
mass definition for m t as we did. 

3) The OPE also shows that the correct scale to be used in the running a s in eq. (33) is the 
high scale Q = 0(m t , M H+ ) and not the low scale m b . 

4) The absorption of the large logarithms into the running mass does not work for terms that are 
suppressed by higher powers (m b /M H+ ) n with respect to the leading contribution considered 
by us. Higher-twist operators contain explicit 6-quark fields. At twist-8 there are operators 
of the form m b H + H~bAb, where A is some Dirac structure. Solving the RG equation for 
these operators yields extra evolution factors in addition to the running mass. These effects 
occur in corrections of order m b /Mf I+ and are certainly only of academic interest. 
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Next consider T(H + tb) for the case m b <C m t <C M H+ : in this limit, another large logarithm, 
logm t /M H+ , appears. Now we have to perform the OPE in two steps. In the first step we again 
match the forward scattering operator to local operators as in eq. (54) at a scale Q\ = 0(M H+ ), 
but we treat the top quark as light, so that the dependence on m t now resides in the matrix element 
( H + | Oi| H + ) rather than in the Wilson coefficient. For simplicity we specify Q\ = M H+ . The 
leading power r° is again represented by the twist-4 operator 0±, yet the corresponding Wilson 
coefficient lacks the factor of (1 — r t ) 2 compared to eq. (55). The terms of order r\ are represented 
by 02 = m\Ox with C 2 = -2C 1 /M 2 J+ . At twist-8 different operators of the form m^H^H^tAt 
with non-trivial anomalous dimensions occur as discussed in point 4 above. In the second step one 
applies an OPE at the scale Q2 = 0{m t ). At this step the dependence on m t migrates from the 
matrix elements into the Wilson coefficients, which at order r\ amounts to a trivial rescaling of the 
coefficients and operators by m t or l/m t . To order r° and r\ the only effect of the OPE is to replace 
the top mass in the expression for T(H + — > tb) in eq. (33) by a running top mass m t (M H+ ), and 
to omit the explicit term proportional to rtlogr< in the 0(a s ) correction. Since we have adopted 
the on-shell definition for the top mass, one must either use a running mass definition based on the 
pole mass (i.e. with mt(m t ) = m? t ole ) or transform the result in eq. (33) to the MS scheme with the 
appropriate change in the 0(a s ) correction. It is a nice check to expand the running mass to first 
order in a s : 

n (M H+ ) = = r t (l + 2 ^ log r t ) + O (afj with r t = r t (m t ) 

and to verify that the overall factor (1 — r t (M H+ )) 2 indeed reproduces the r t logr t term in eq. (33). 
The terms of order r 2 log r t are not correctly reproduced by the running top mass as anticipated 
by the occurrence of non-trivial twist-8 operators. The important result of our consideration of the 
case m b <C m t <C M H+ is the absence of terms of the form r° log r t to all orders in a s . In this case 
the additional large logarithm log r t is always suppressed by powers of r t and therefore these terms 
are negligible for r t <C 1 and need not be resummed. 

For the decay t — > bH + the above discussion can be repeated with the appropriate changes in 
the OPE: the leading-twist operator is now 0\{Q) = m 2 {Q)tt{Q) and the external state in eq. (54) 
is a top quark instead of a charged Higgs boson. We have m b <C m t ,M H+ and the factorization 
scale Q is again of order m t ,M H+ . While 0\ now involves strongly interacting fields, its matrix 
element {t \0± \ t)(Q) still does not contain large logarithms log m b /Q other than those contained in 
the running mass m^Q). Hence the proof above for T{H + — ► tb) applies likewise for T(t — ► bH + ). 

After exchanging Vt^tibR for V&ChbR in eq. (21), we can likewise apply our effective lagrangian C 
to semileptonic i?-meson decays corresponding to b — > ctT>i by using the appropriate scale Q ~ m b 
in C. The QCD radiative corrections involve no large logarithm, because the gluons couple only to 
the b and c quarks. Hence the effective four-fermion operator ci,bji ~t-RVL obtained after integrating 
out the heavy H + renormalizes in the same way as the quark current c_l&_r in C. The corresponding 
loop integrals do not depend on M R+ at all and this feature is correctly reproduced by using mb(m b ) 
in C The situation is different in physical processes in which the charged Higgs connects two quark 
lines, as for example in the loop-induced decay 6^57. Here effective four-quark operators, which 
involve a non-trivial renormalization group evolution, occur. The large-tan/3 supersymmetric QCD 
corrections associated with Am b and the H + tLbR Yukawa coupling, however, are still correctly 
reproduced by applying C to b — > s 7 or other loop-induced rare 6-decays. Yet it must be clear that 
these corrections are part of the mixed electroweak-QCD two-loop contributions and that there are 
already supersymmetric electroweak contributions at the one-loop level, which are process-specific 
and of course not contained in C. 
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